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We evaluated genetic variation in Middle East respiratory 
syndrome coronavirus (MERS-CoV) imported to South Ko- 
rea in 2018 using specimens from a patient and isolates 
from infected Caco-2 cells. The MERS-CoV strain in this 
study was genetically similar to a strain isolated in Riyadh, 
Saudi Arabia, in 2017. 


etween 2012 and the end of August 2018, a total of 

2,248 laboratory-confirmed cases of Middle East re- 
spiratory syndrome coronavirus (MERS-CoV) infection 
and 798 associated deaths (case-fatality rate 35.3%) were 
reported in 27 countries (/). Most of these cases (83%) 
were reported in Saudi Arabia (1,871 cases, including 724 
related deaths, for a case-fatality rate of 38.7%). 

The 2015 MERS-CoV outbreak in South Korea re- 
sulted in 186 laboratory-confirmed cases and 38 deaths. 
MERS-CoV spread via intrahospital and interhospital 
transmission at 16 clinics and hospitals across South Korea 
(2). This outbreak was recorded as the largest MERS out- 
break outside of the Middle East (3). 

In August 2018, eight laboratory-confirmed cases of 
MERS-CoV infection were reported in Saudi Arabia, with 
4 associated deaths; | case in a person who had travel his- 
tory to the Arabian Peninsula was reported in the United 
Kingdom (4). In September 2018, a man in South Korea 
with a history of travel to the Middle East became ill and 
was suspected of having MERS-CoV infection. We report 
the investigation of this case and genetic characterization 
of the virus on the basis of variation and phylogenetic anal- 
yses of the spike (S) gene and the full-length genome from 
patient specimens and isolates from Caco-2 cells. 
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The Study 

On September 7, 2018, a 61-year-old man who had moved 
from Kuwait to South Korea visited the hospital with diar- 
thea and fever (5). The patient had a history of visiting the 
Middle East and was suspected of MERS-CoV infection; 
accordingly, he was transferred to a hospital with a national 
isolation ward. 

On September 8, we collected nasopharyngeal swab 
and sputum samples from the patient and transported 
them to the Seoul Provincial Institute of Public Health 
and Environment (Seoul, South Korea). Real-time re- 
verse transcription PCR showed that the sputum was 
positive for the upstream regions of the E protein (upE) 
and orfla genes (6). Later the same day, we collected 
nasopharyngeal swab, throat swab, and sputum samples 
from the patient for further confirmation. The Korea 
Centers for Disease Control and Prevention (KCDC) 
confirmed that throat swabs and sputum samples were 
positive for upE and orfla. 

We added throat swab and sputum samples from the 
patient to a monolayer of Caco-2, Huh 7, and Vero E6 cells, 
which we incubated at 37°C in 5% CO,. We observed a 
cytopathic effect on the third day after inoculation in Caco- 
2 cells and confirmed viral replication of MERS-CoV in 
the supernatant of the infected cells by real-time reverse 
transcription PCR. 

We extracted viral RNA from the specimens and the 
supernatant of the Caco-2 cells using the Viral RNA Mini 
Kit (QIAGEN, http://www.qiagen.com). We performed 
cDNA synthesis using the Superscript IV First-Strand 
synthesis system (Thermo Fisher Scientific, http://www. 
thermofisher.com) with random hexamers. We ampli- 
fied cDNA by overlapping PCR to generate products of 
600—1,100 bp covering the entire S gene and the full-length 
genome; we sequenced the resulting PCR amplicons by 
Sanger sequencing using an ABI 3730 Analyzer (Applied 
Biosystems, http://www.thermofisher.com). 

We used ClustalW  (http:/(www.clustal.org) to 
align the S gene sequences of the isolates from the pa- 
tient, MERS-CoV/KOR/KCDC/001_2018/SP-1 (from 
sputum sample), TS-1 (from throat swab sample), and 
TSVi (from Caco-2 processing), with those of 157 oth- 
er MERS-CoV strains from GenBank and the MERS- 
CoV sequence database (https://www.ncbi.nlm.nih.gov/ 
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Figure 1. A portion of the molecular phylogenetic tree and coding region variants for the spike (S) glycoprotein gene of MERS-CoV 
isolates from South Korea, September 2018, and reference sequences. Phylogenetic analysis of 157 S sequences was performed using 
MEGA7 (https:/Awww.megasoftware.net), with tree visualization using FigTree version 1.4.3 (http://tree.bio.ed.ac.uk/software/figtree). 
The taxonomic positions of circulating strains from the outbreak in South Korea and Riyadh, Saudi Arabia, are indicated. Boldface 
indicates compressed major clades of MERS-CoV. Bootstrap values (>70%) on nodes are shown as percentages based on 1,000 
replicates. Red indicates South Korea 2018 isolates, blue indicates South Korea 2015 isolates, and green indicates Saudi Arabia 2017 
isolates. The highlighted strain, from 2012, is the prototype strain used as the reference. EMC, Erasmus Medical Center; MERS-CoV, 
Middle East respiratory syndrome coronavirus. Scale bar indicates nucleotide substitutions per site. An expanded version of this figure 
showing the full phylogenetic tree is available online (https:/Awwwnc.cdc.gov/EI|D/article/25/5/18-1534-F 1.htm). 


genomes/VirusVariation/Database/nph-select.cgi?cmd= 
database&ctaxid=1335626), deposited up to August 2018. 
We sequenced and analyzed the full-length genome of 
MERS-CoV/KOR/KCDC/001_2018/TSVi and 39 other 
MERS-CoV strains. We constructed phylogenetic trees by 
the maximum-likelihood method with 1,000 bootstrap rep- 
licates using MEGA7 (7) and RAxML (https://github.com/ 
stamatak/standard-RAxML). We constructed phylogenetic 
trees from the MERS-CoV S genes (4,062 bp) (Figure 1) 
and the full-length genome (30,150 bp) of the virus (Figure 


Emerging Infectious Diseases * www.cdc.gov/eid * Vol. 25, No. 5, May 2019 


2) obtained from the patient in this study and the most simi- 
lar human MERS-CoV sequences from other countries (4). 

We submitted viral S gene sequences MERS-CoV/ 
KOR/KCDC/001_ 2018 SP-1, TS-1, TSVi from the patient 
to GenBank under accession nos. MH978886, MH978887, 
and MH978888. MERS-CoV/KOR/KCDC/001_2018 
SP-1, TS, and TSVi were most closely related (99.85%— 
99.90% nucleotide identity) to a human MERS-CoV strain 
isolated in Riyadh, Saudi Arabia, in 2017 (Hu-Riyadh- 
KSA-9730_2017; GenBank accession no. MG912608). 
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KOR/KCDC/2018 SP-1 and TS had 5 nucleotide substi- 
tutions (C309T, TS519C, C2928T, T3375C, and T3598C) 
compared with a MERS-CoV strain found in Riyadh in 
2017, whereas MERS-CoV/KOR/KCDC/001_ 2018 TSVi 
had 3 nucleotide substitutions (C309T, C2928T, and 
T3375C) compared with the same strain. There were no 
substitutions in the N terminal domain (NTD) and recep- 
tor-binding domain (RBD) of the S gene in the same com- 
parison (Table). We submitted the full-length genome se- 
quences from MERS-CoV/KOR/KCDC/001_2018 TSVi, 


obtained from cell isolates, to GenBank under accession 
no. MK129253. In a phylogenetic analysis, the virus was 
located in a lineage of locally endemic 2016 and 2017 
Saudi Arabia strains (lineage 5 in Figure 2) but was distin- 
guished from 2015 South Korea strains. 


Conclusions 

We characterized the S gene and full-length genome of 
MERS-CoV imported to South Korea from the Mid- 
dle East using sputum and throat swab specimens from 


H310912 | Saudi Arabia | Hu/Riyadh-KSA-7344/2017 | 2017/05/15 
MG011352 | Saudi Arabia | Hu/Qasim-KSA-893/2016 | 2016/07/20 
MG011359 | Saudi Arabia | Hu/Riyadh-KSA-16559/2016 | 2016/03/29 
MH395139| Saudi Arabia | 2363 |2016/12/04 

MG912608 | Saudi Arabia | Hu/Riyadh-KSA-9730/2017 | 2017/06/14 
MG912607 | Saudi Arabia | Hu/Riyadh-KSA-9522/2017 | 2017/06/11 
MG011361| Saudi Arabia | Hu/Riyadh-KSA-K17000405/2017 | 2017/01/10 
MK129253| Korea| MERS-CoV/KOR/KCDC/001-2018-TSVi| 2018/10/20 
MG011354 | Saudi Arabia | Hu/Riyadh-KSA-20570/2016 | 2016/06/16 

KX154692 | Saudi Arabia | Hu/Riyadh-KSA-12969/2016 | 2016/02/27 

KX154694| Saudi Arabia | Hu/Artawiyah-KSA-13328/2016| 2016/02/29 
MH306207 | Saudi Arabia | HCoV-EMC|2016/08/01| No clade2 

MG011342 | Saudi Arabia | Hu/Riyadh-KSA-13984/2016| 2016/03/07 
MG011351| Saudi Arabia | Hu/Madinah-KSA-390/2016 | 2016/04/10 
KT225476|Thailand| MERS-CoV/THA/CU/17 06 2015|2015/06/17 

KT806054| Saudi Arabia | Hu/Najran-KSA-C20915/2015 | 2015/02/13 

KT806044| Saudi Arabia | Hu/Jeddah-KSA-C20843/2015 | 2015/02/09 

88" KT806045 | Saudi Arabia | Hu/Jeddah-KSA-C21271/2015 | 2015/02/22 

KT036372| China| ChinaGDO1|2015/05/28 

KT326819 | Korea| MERS-CoV/KOR/KNIH/001 05 2015 | 2015/05/22 

98 '— KT029139 | Korea| MERS-CoV/KOR/KNIH/002 05 2015 | 2015/05/20 


2 KF600628| Saudi Arabia |Hafr-Al-Batin 1 2013 |2013/06/04 
KT368825 | Saudi Arabia | Camel/Riyadh/Ry23N/2014 | 2014/07 
KM027262 | Saudi Arabia | Riyadh 2014KSA 683/KSA/2014| 2014/04/22 


KM027257|Saudi Arabia |Jeddah C7770/KSA/2014-04-07 | 2014/04/07 

KM027260|Saudi Arabia | Jeddah C10306/KSA/2014-04-20|2014/04/21 

KRO11263| Saudi Arabia | Hu/Riyadh-KSA-2345/2015 | 2015/01/21 

KJ713295 | Saudi Arabia | KSA-CAMEL-505 | 2013/11 

100] KM015348| United Kingdom | England/2/2013|2013/02/10 

KM210278| United Kingdom | England/3/2013|2013/02/10 

KF186564| Saudi Arabia|Al-Hasa 4 2013| 2013/05/01 

KP209306| United Arab Emirates | Abu Dhabi UAE 8 2014|2014/04/07 

KP719933| United Arab Emirates | Camel/UAE/D1209/2014 | 2014 

100, KF917527|Saudi Arabia| MERS-CoV-Jeddah-Camel-1| 2013/11/08 
KF958702|Saudi Arabia | MERS-CoV-Jeddah-human-1| 2013/11/05 

KJ156949| Saudi Arabia|Taif 1 2013|2013/06/12 

KJ156881 | Saudi Arabia| Wadi-Ad-Dawasir 1 2013 |2013/06/12 

KC776174| Jordan |Jordan-N3/2012|2012/04 

100 JX869059 | | Human betacoronavirus 2c EMC/2012|2012/06/13 


100 


100 


Lineage 3 


Lineage 4 


Lineage 2 


100 


Lineage 1 


0.0005 
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Figure 2. Maximum-likelinood tree showing MERS-CoV isolates from South Korea, September 2018 (black dot), and reference 
MERS-CoV genomes. Tree estimated using RAxML values (https://github.com/stamatak/standard-RAxML) on branches are shown 
as percentages based on 1,000 bootstrap replicates from the nucleotide sequences. MERS-CoV, Middle East respiratory syndrome 


coronavirus. Scale bar indicates nucleotide substitutions per site. 
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Table. Comparison of spike glycoprotein gene sequence variants for a MERS-CoV isolates from a patient in South Korea, 2018, and 


reference strains* 


Nucleotide (amino acid) positions 


NTD RBD Other regions 

183 258 409 1586 1588 1848 1886 3177 3267 
MERS-CoV strain (61) (86) (137) (529) (530) (616) (629) (1059) (1089) 
JX869059_HCoV-EMC_2012t C(G) T(V) A(S) T(l) G(V) T(VY) G(R) T(D) C(S) 
KRO11266_KSA_Hu/Riyadh-KSA-2049/2015 C (Vv) 
KRO11263_KSA_Hu/Riyadh-KSA-2345/2015 C (Vv) 
KRO11264_KSA_Hu/Riyadh-KSA-2343/2015 C (Vv) 
KR0O11265_KSA_Hu/Riyadh-KSA-2466/2015 C (Vv) 
KT026453_KSA_Hu/Riyadh-KSA-2959/2015 C (Vv) 
KT026454_KSA_Hu/Riyadh_KSA_4050/2015 C (Vv) 
KX154690_Hu_Jeddah-KSA-161RS1146_2016 C (Vv) 
KX154684_Hu-Riyadh-KSA-11739_2016 C (Vv) 
KX154694_Hu-Artawiyah-KSA-13328_2016 C (Vv) 
MG011361_Hu-Riyadh-KSA-K17000405_2017 C (Vv) 
MG011360_Hu-Riyadh-KSA-K17000887_2017 C (Vv) 
MG912608_Hu-Riyadh-KSA-9730_2017 C (Vv) 
MERS-CoV/KOR/KNIH/001_05_2015 C (Vv) C(T) C(D) T(S) 
KT029139-MERS-CoV/KOR/KNIH/002_05_2015 cC(Vv) C(R) C (L) C(D) T(S) 
MERS-CoV/KOR/KNIH/009_05_2015 C (Vv) C (Vv) C(D) T(S) 
KT006149-ChinaGD01_ 2015 G(G) C(V) C (Vv) C(D) T(S) 
MERS-CoV/KOR/KNIH/012_05_2015 C (Vv) C (T) C (Vv) C(D) T(S) 
MERS-CoV/KOR/KNIH/013_05_2015 C (Vv) C (T) C (Vv) C(D) T(S) 
MERS-CoV/KOR/KNIH/015_05_2015 C (Vv) C(T) C (Vv) A(H) C(D) T(S) 
MERS-CoV/KOR/KNIH/042_04_2015 C (Vv) C(D) T(S) 
MERS-CoV/KOR/KCDC/001_2018_2SP-1 C (Vv) 
MERS-CoV/KOR/KCDC/001_2018_4TS C (Vv) 
MERS-CoV/KOR/KCDC/001_ 2018 _4TSVi C (Vv) 


*Only variations from the reference strain are shown. Boldface indicates the South Korea isolates and the most similar strain identified in this study. MERS- 
CoV, Middle East respiratory syndrome coronavirus; NTD, N-terminal domain; RBD, receptor-binding domain. 


tReference MERS-CoV strain. 


a patient and from Caco-2 cells exposed to throat swab 
samples from the patient (89). A comparative genet- 
ic analysis of 157 strains isolated in 27 countries in- 
cluding South Korea demonstrated that the imported 
MERS-CoV we isolated showed the highest identity 
(99.85%-99.90%) with a strain recently isolated in Saudi 
Arabia in 2017 (/,/0). These strains differed at 2 positions 
(C309T and T519C). We detected no variation in NTD or 
RBD, which are the major functional sites of the S protein 
(71). In addition, the strains involved in the outbreak in South 
Korea in 2015 had 12 nt substitutions (NTD, | site; RBD, 2 
or 3 sites), implying that the imported case in 2018 involved 
a different virus (4). 

Despite 110 laboratory-confirmed cases of MERS- 
CoV reported globally during 2018, no genetic information 
has been obtained to date. Therefore, we could not compare 
the genetic relationships between the strain we isolated and 
MERS-CoV strains in the Middle East (/2). 

In summary, we obtained genetic information about a 
MERS-CoV strain in South Korea in 2018 that was prob- 
ably imported from the Middle East. We expect the se- 
quence data we identified to support future epidemiologic 
investigations of MERS-CoV, particularly to fill the need 
for recent data. Furthermore, the results of this study im- 
prove our understanding of the evolution of MERS-CoV 
and highlight the need for enhanced surveillance, espe- 
cially for persons traveling to the Middle East. 
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